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We discuss the QCD phase diagram from two different point of view. 
We first investigate the phase diagram structure in the strong coupling 
lattice QCD with Polyakov loop effects, and show that the the chiral and 
Zn c deconfincment transition boundaries deviate at finite y, as suggested 
from large N c arguments. Next we discuss the possibility to probe the 
QCD critical point during prompt black hole formation processes. The 
thermodynamical evolution during the black hole formation would result in 
quark matter formation, and the critical point in isospin asymmetric matter 
may be swept. {T^b) region probed in heavy-ion collisions and the black 
hole formation processes covers most of the critical point locations predicted 
in recent lattice Monte-Carlo simulations and chiral effective models. 

PACS numbers: 26.50. +x, 25.75.Nq, 12.38.Gc, ll.10.Wx, 11.15.Me, 12.39.Fe 

1. Introduction 

QCD phase transition at finite temperature (T) and chemical potential 
(fj,) is attracting much attention in recent years. The beam energy and 
system size scan programs at RHIC pQ and SPS [2] are running to discover 
the critical point and the first order transition at finite T and fj,, and the 
discovery of the two solar mass neutron star [3] gives us some hints on the 
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phase transition at finite density. Recent large N c arguments suggest the 
existence of another form of matter, referred to as the quarkyonic matter, 
where the Polyakov loop is suppressed and the density is high [3]. 

In this proceedings, we first discuss the chiral and Zn c deconfinement 
transitions at finite T and fi based on the strong coupling lattice QCD with 
Polyakov loop effects (P-SC-LQCD) [5]. Next we discuss the possibility to 
probe the critical point in isospin asymmetric high baryon density matter 
formed during the prompt black hole formation processes [6]. 

2. Strong coupling lattice QCD and quarkyonic matter 

Do the chiral and Z^ c deconfinement transition boundaries deviate at 
large fi ? This is one of the most interesting questions in the current QCD 
phase diagram studies. From the large N c arguments, we expect the exis- 
tence of the so-called quarkyonic matter, where the Zn c order parameter 
(Polyakov loop) is suppressed and the density is high [3J. In the lattice 
QCD Monte-Carlo simulations at fi = 0, the Zn c transition temperature 
(T c (Zn c )) is close to but somewhat larger than the chiral transition tem- 
perature (T c (x)) [3- I n chiral effective models, some of them predict the 
existence of quarkyonic-like matter [SJ [H] , while some of them predict that 
the Zn c transition boundary agrees with the chiral transition boundary |lUj . 
Thus it is important to discuss the chiral and Z^ c transition boundaries at 
finite [i in the theoretical framework directly based on QCD for N c = 3. 

The strong coupling (1/g 2 ) expansion in the lattice QCD (SC-LQCD) 
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Fig. 1. Left: Comparison of T c at fi = in P-SC-LQCD (solid) and SC-LQCD 
(without the Polyakov loop, dashed) [5]. Two different Polyakov loop treatments 
(Weiss and Haar method) are compared. The triangles represent the lattice MC 
results with one species of unrooted staggered fcrmion. Right: The phase boundary 
for the chiral transition with the second peak of d£ p /dT at j3 = 4 in the chiral limit. 
The "P" and "Q" correspond to those in Fig. [2J 
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has been successful since the beginning of the lattice gauge theory [11] and 
would provide an alternative lattice framework to study the QCD phase 
diagram at finite T and \x [121 113] . For example, the spontaneous breaking 
of the chiral symmetry and its restoration at finite T and/or [i have been 
known to be realized in the strong coupling limit [13| . Recently, the finite 
couping and Polyakov loop effects are incorporated in the framework of SC- 
LQCD, and are found to explain the MC results of T c at \x = in the region 
P = 2N c /g 2 < 4 13 HI|, as shown in the left panel of Fig. [TJ 

We discuss here the chiral and Z^ c deconfinement dynamics by using the 
SC-LQCD with the Polyakov loop effects (P-SC-LQCD) in the mean field 
approximation [5]. We take account of the next-to- leading order (NLO, 
0(l/g 2 )) and the leading order (0(1/ g 2Nr )) of the strong coupling expan- 
sion in the fermionic and pure gluonic sector, respectively, and in the leading 
order of the 1/d expansion [15]. The effective potential is given as [5] 

J- Cff ($; T, M ) EE -(TlogZ LQ CD)/< = + -Feff 01 , (1) 

^ - + f>.<P.) ^ + ^ + y (<P 2 r ~ <4) ~ N c bg ^Z~ZZ 

- N c E q - T [log TZ(E q - fi, £, I) + log lZ(E q + ft, I, £)] , (2) 
jf ff o1 ~ -2TdN 2 c (-^f) " Tlog Mn^ilp, i P ) , (3) 

jX = fj, - log ^Z + /Z_ , Z± = 1 + P T (ip T ± UJ T ) , 

1l(x, £, l) = l + N c ie~ x l T + N c £e- 2x / T + e~ Zx ' T , (4) 
Ma*(f, l) = l- W - m? + 4(^ 3 + P) , (5) 

where a and £(£) denote the chiral condensate and the (anti-)Polyakov 
loop, respectively, d = 3 is the spatial dimension, f3 T = /3d/2N^ and 
(3 S = pd(d — 1)/16A^. Several other auxiliary fields (ip s , ip T ,Lu T ) are in- 
troduced in addition to a and £, and the stationary conditions are imposed 
on these fields in equilibrium. The Polyakov loop £ couples with quarks, and 
appears with the Boltzmann factor e~^ Eq ~^^ T . Color-singlet states domi- 
nate in the confined phase (I ~ 0), while quarks can excite in the deconfined 
phase (I 0). This point has been pointed out in the strong coupling limit 
and utilized in the PNJL model |16j . 

We discuss the results at P = 4, which is in the coupling range where 
P-SC-LQCD roughly explains the LQCD-MC results of T c at fi = 0, and 
results are shown in the lattice unit, a = 1. In the right panel of Fig. [TJ 
we show the phase diagram at /3 = 4 in the chiral limit (mo = 0). We 
find the second- (dashed) and first-order (solid) chiral transition boundaries 
separated by the (tri-)critical point (CP) at (/UcPj^cp) — (0.58,0.19). 
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Fig. 2. Left: T dependence of the chiral condensate (dashed) and the Polyakov 
loop (solid) at /i = 0.5 and /3 = 4 in the chiral limit. "P" and "Q" correspond to 
the points in the right panel of Fig. [T] Right: T dependence of da/dT (dashed) 
and d£ p /dT (solid) with m Q = 0.03 at /3 = 4 and fi = 0.5. 



In the left panel of Fig. [2j we show the T dependence of the chiral 
condensate a and the Polyakov loop £ p at /U = 0.5. At the chiral transition 
temperature, the chiral condensate becomes zero quickly, and the Polyakov 
loop is affected to have a kink. In addition to this chiral-induced kink shown 
as "P" in the figure, we find another temperature shown as "Q" where the 
Polyakov loop increases rapidly. This feature is more clearly seen in d£ p /dT. 
We show d£p/dT at /i = 0.5 and m-o = 0.03 in the right panel of Fig. [2j We 
can see the second peak in d£ p /dT at Tq ~ 0.52. A similar double-peak 
structure has been reported in the model studies based on PNJL model [9]. 

The peak "Q" can be understood as the Z/v c -induced crossover from 
following reasons. The temperature Tq ~ 0.52 is found to be almost inde- 
pendent on the chemical potential, as indicated by the upper line in the right 
panel of Fig. [TJ Tq is also insensitive to quark mass. The temperature of the 
peak "P" is shifted upward and becomes closer to "Q" with increasing mo, 
while Tq stays almost constant. For larger masses, mo > 0.05, the two peaks 
merges to one, and it grows with increasing mo at nearly mo independent 
temperature. These observations agrees with the expected character of the 
Z7v c -induced transition; the Z^ c deconfmement transition nature would be 
stronger with large quark mass, and its transition temperature would have 
small dependence on /i from the large N c argument. This interpretation 
supports the existence of the quarkyonic-like matter in cold-dense matter, 
where the Polyakov loop is suppressed and density is high. 

3. Critical Point Sweep during Black Hole Formation 

Critical point (CP) is one of the largest targets in the beam energy and 
system size scan programs at RHIC and SPS and in the forthcoming FAIR 
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Fig. 3. The BH formation profile, (T, ^j,b, as a function of the radius. Results 
are shown for the gravitational collapse of a 40 Mq star at t = 0.5 sec (dotted), 
1.0 sec (dashed), and 1.344 sec (solid lines, just before the horizon formation). 



facility. CP is expected to be probed in these experiments if it is in the low 
baryon chemical potential region, /iqp < 500 MeV, as predicted in some of 
the lattice MC calculations [26 1 127 1 128]. On the other hand, we cannot reject 
the possibility that the CP is located in the lower T and higher jjl region, 
as predicted in many of the chiral effective models [2TJ CE21 E21 EH IM]- 
Therefore, it is important to search for other candidate sites where hot and 
dense matter is formed and CP is reachable. 

A gravitational collapse of a massive star is a promising candidate. A 
majority of non-rotating massive stars with mass M > 2OM are expected 
to collapse quietly (faint-supernova) to black holes (BH) [18] . Their fre- 
quency should be comparable to supernovae provided that the mass spec- 
trum of stars has a long tail as in the power law behavior. The BH formation 
processes are found to form hot (T ~ 90 MeV) and dense {ps ~ 4po) mat- 
ter in the neutrino-radiation hydro dynamical simulations in the collapse 
and bounce stage of a 40-M© star [19]. Thermodynamical variables at a 
given time vary as a function of radius in a proto-neutron star and form a 
line (referred to as the BH formation profile in the later discussions) in the 
T—pB plane. In Fig. [3l we show the BH formation profile (T, /ie,<5/i) [19] 
calculated by using the Shen EOS at t = 0.5,1.0 and 1.344 sec after the 
bounce during the BH formation from a 40 Mq star in the proto-neutron 
star core, where the mass coordinate from the center is M < 1.6M & . The 
time t = 1.344 sec is just before the horizon formation. From the outer to the 
inner region of the proto-neutron star, T first increases from T ~ 10 MeV 
to T ~ (50 — 90) MeV in the middle heated region, and decreases again 
inside. The baryon chemical potential hb is found to go over 1300 MeV 
in the central region just before the horizon formation at t = 1.344 sec. 
The isospin chemical potential is found to be <5/i = (50 — 130) MeV in the 
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inner region. The temperature and density in BH formation are signifi- 
cantly higher than in the model explosion calculation of supernova. The 
highest temperature and density are moderate in the collapse and bounce 
stage of supernovae, (T,ps) ~ (21.5 MeV,0.24 fm~ 3 ) when hadronic EOS 
is adopted |20] , while it has been argued that the transition to quark matter 
might trigger successful supernovae [T?] . 

We discuss here the possibility that the BH formation profile evolves 
with time and may pass through CP and the vicinity (CP sweep). The 
CP location scatters in the T—jjLB plane in chiral effective models such as 
NJL [21], P-NJL [EJ [22], P-NJL with 8-quark interaction (P-NJL 8 ) [23], 
and PQM [21] models. We expect that CP moves in the lower T direction 
at finite 5p, because d-quark dominates and the effective number of flavors 
decreases. Since the matter passes through the high ps an d low T region 
compared with high-energy heavy-ion collisions, the reduction of the CP 
temperature Tcp is essential for the CP sweep during the BH formation. 

The Lagrangian density of the Polyakov loop extended quark meson 
(PQM) model, as an example of chiral effective models, is given by 

£ =q [ij^D^ - g(o + ij 5 T ■ tt) - g v ^{u^ + r • R^)] q 

+\{d^f + \{d^f ~ U(a, tt) - Uiil, I, T) 

- ^JV • + lrn 2 v (uj^ + R^ ■ RP) . (6) 

The mesonic potential is U(a, 7r) = A (<r 2 + 7v 2 — u 2 ) 2 /4 — ha, and oj^ and 
R^v are the field tensors of the uj and p mesons. We use the Polyakov loop 
effective potential Ue(£,£,T) = T 4 [-a(T)U/2 + b(T) log H{1 ,£)] where the 
Polyakov loop is defined as t = Ti[Pexp(i J^' T dTA/^)\/N c . In the mean 
field approximation, the thermodynamic potential is obtained as 




where E p = \/ p 2 + M 2 , M = ga + mo is the constituent quark mass, and 
jlf = p =F SfM — 2g 2 pj /m 2 is the effective chemical potential with qp = — , + 
for f = u and d, respectively. While the PQM model is renormalizable, we 
adopt here the momentum cutoff A for simplicity. 

We show the first order phase boundary of symmetric (5p = 0) and 
asymmetric (5p ^ 0) matter in PQM in the left panel of Fig. HJ We find a 
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Fig. 4. Left: First order phase boundary surface (solid lines) in (T, hb,Sii) space 
calculated with the PQM model are compared with the BH formation profile, (thcr- 
modynamical profile (T, hb, 8fi) during the BH formation) at t = 0.5, 1.0, 1.344 sec 
after the bounce (double lines). Right: Predictions of the critical temperature and 
critical point locations in comparison with the chemical freeze-out points and the 
swept region during the black hole formation. 



trend that the first order phase boundary shrinks at finite 5ji. Transition 
temperature at a given baryon chemical potential fiB = decreases, and 
the transition baryon chemical potential /i c at T = also decreases. We 
do not consider here the pion condensed phase, because the s-wave ttN 
repulsion would suppress the s-wave pion condensation [25]. 

The CP location is sensitive to <5/z. Compared with the results in 
symmetric matter, Tcp becomes smaller at finite 5fi and reaches zero at 
5/j, = 5fi c ~ (50 — 80) MeV. CP is also sensitive to the model and param- 
eters as shown in the right panel of Fig. [H The Polyakov loop suppresses 
single quark excitations in the hadron phase, then the transition temper- 
ature and thus the critical temperature are shifted upwards in the P-NJL 
model. The temporal component of the vector potential shifts the chemical 
potential effectively, and consequently leads to an upward shift of /i c by 
about 10-15 MeV at g v /g = 0.2. 

We shall now compare the CP location and the phase boundary with 
the BH formation profile. In the left panel of Fig. 01 we compare the phase 
boundaries and the BH formation profile in the PQM model. During the 
BH formation, the baryon chemical potential reaches around 1000, 1100 
and 1300 MeV in the central region of the proto-neutron star at t = 0.5, 1.0 
and 1.344 sec, respectively, suggesting that quark matter would be formed 
between t = 0.5 and 1.0 sec in the central region in most of the models 
considered here. We also find that the BH profile go through the critical 
line in asymmetric matter, i.e. the CP sweep takes place in PQM. In other 
models, we find that there are three possible types in the transition to the 
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quark matter during the evolution of matter toward the BH formation; the 
first order transition, the crossover transition, and the CP sweep, where 
the BH formation profile goes below, above and through the critical line in 
asymmetric matter. 

The predicted CP locations in lattice MC [HJEIIEE] anc i effective mod- 
els [6] seem to be in the (T,/j,b) region probed in heavy-ion collisions [291130] 
or during the prompt black hole formation processes [6] as shown in the right 
panel of Fig. [5J Recent LQCD-MC predictions by using the reweighting 
(LR04) [26], Taylor expansion (LT04) [27], and canonical ensemble method 
(LC11) |28] are consistent with each other, and suggest that the CP is 
accessible in heavy- ion collisions. It should be noted that some previous 
studies [3TJ predicted larger ^cp> and there is also implication that there is 
no chiral critical point in the small fi region for Nf = 3 [32] . By comparison, 
effective models generally predict the CP in the lower T and larger fi region, 
and many of them are accessible during the black hole formation. Further 
studies are necessary to understand the difference between the lattice MC 
and effective model results. 

4. Summary and Discussion 

We have discussed here two subjects on the QCD phase diagram at finite 
density. One of them is the chiral and Z^ c deconfinement transition bound- 
aries at finite fx. In the strong coupling lattice QCD with Polyakov loop 
effects [5], the Z^ c deconfinement boundary defined as the peak of di/dT 
is found to deviate from the chiral transition boundary at finite /x, and it 
suggests the existence of the Polyakov loop suppressed high density matter, 
which may be interpreted as quarkyonic matter [3]. This would be consis- 
tent with the lattice QCD Monte-Carlo simulation results [7J, which suggests 
that the Zn c deconfinement transition temperature T c (Pol.) would be mean- 
ingfully higher than the chiral transition temperature T c (x) at fi = 0. The 
chemical potential effects should be stronger on T c {x) than on T c (Pol.), then 
separation at fj, = (T c (x) > T c (Pol.)) would be enhanced at finite fi. This 
contradicts to the results including additional chiral-Polyakov coupling, or 
the chemical potential dependence of the Polyakov loop potential [ID]. In 
the strong coupling lattice QCD, the additional chiral-Polyakov coupling 
appears in the higher order terms of the large-dimensional (1/d) expansion, 
and its effects should be studied. 

The possibility to probe the QCD critical point during black hole forma- 
tion processes is discussed in the second part [6j. We have found that the 
critical point in isospin asymmetric matter would be accessible in the black 
hole formation processes, if the CP is in the low T and high /x region as 
predicted in chiral effective models. It should be noted that we compare the 
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results of the CP location in chiral effective models and the thermodynam- 
ical condition (T, /xb) calculated with the hadronic EOS. This comparison 
is relevant, since the thermal trajectory should be the same even if we use 
the combined EOS of quark and hadronic matter, as long as the hadronic 
EOS is reproduced at low T and in the combined EOS. It is desired 
to examine the thermodynamical evolution in the combined EOS. Another 
interesting point is the differences of the CP locations in chiral effective 
models and in lattice MC simulations. While lattice MC simulations have 
the sign problem and there is also implications that there is no chiral criti- 
cal point in the low fi region [32] , recent results consistently suggest the CP 
may exist in the low chemical potential region, ^cp — 500 MeV. It would 
be necessary to understand these differences in order to pin down the CP 
location in the phase diagram. 
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